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5-H ^ This paper investigates the evolution of a void and an adjacent galaxy supercluster. 

'"H, For this purpose the quasispherical Szekeres model is employed. The Szekeres model 

■^r ' is an exact solution of the Einstein field equations. In this way investigations of the 

evolution of the cosmic structures presented here can be freed from such assumptions 
as small amplitude of the density contrast. Studying the evolution of a void and an 
adjacent supercluster results not only in better understanding the evolution of cosmic 
structures but also in acquainting us with the Szekeres model. The main results include 
the conclusion that small voids surrounded by large overdensities evolve slower than 
large, isolated voids do. On the other hand, large voids enhance the evolution of adjacent 
supcrclusters which evolve much faster than isolated galaxy superclusters. 

o 
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i ' 1. Introduction 

The structures which can be observed in the local Universe include small voids 
^ ' among compact clusters, superclusters and large voids surrounded by large walls 

OO . or long filaments. The present day density contrast of overdense regions is larger 

J^^ \ than 1^ and inside voids it descends to -1.^ To describe these structures, an exact 

solution of the Einstein field equations must be employed. However, among the 

known exact solutions none is flexible enough to describe such complicated structure 
f^ ' as our Universe. Therefore, the analysis of this paper will focus on smaller scales. 

52 ' The structures on small scales, up to Mpc can be described by the Szekeres model 

which is an exact solution of the Einstein field equations. 

The structure of this paper is as follows: Sec. [2] presents the Szekeres model; Sec. 

[3] presents the evolution of pairs void-supercluster in the quasispherical Szekeres 
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5h ' model; Sec. |4] presents the role of expansion in the process of structure formation. 

2. The Szekeres model 

For our purpose it is convenient to use a coordinate system different from that in 
which Szekeres^ originally found his solution. The metric is of the following form:^ 

as cat ^^_ ^^ ar ^ ^2 ' i-^i 

where ' = d/dr, e — ±1, and k = k{r) < e is an arbitrary function of r. 
The function E is given by: 
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where the functions S = S{r), P — P{r), Q = Q{r), and C — C{r) satisfy the 
relation: 

p2 (Q2 c 

^■=25+l + r- = = »•"> W 

but are otherwise arbitrary. 

The quasispherical Szekeres model is the case of e = 1. As was shown in Ref. 5, 
surface of t = const and r — const has a topology of a sphere. However, as S, P 
and Q are now functions of r, the spheres are not concentric. For the spheres to 
be concentric, the functions S, P and Q must be constant. Such conditions entail 
spherical symmetry, with which the Szekeres model becomes the the Lemaitre- 
Tolman model. ^'^° 

Substituting the metric ([T]) in the Einstein equations, and assuming the energy 
momentum tensor for a dust, the Einstein equations reduce to the following two: 

72^'(^'^^^='^)-Hr) + lA<^Ht,r), (4) 



4n^o(t rva)^ M'{r) - 3M{r)E'{r,p,q)/E{r,p,q) 

^2/^V ' "^'^^ ^^{t,r)[<^'{t,r)-<i>{t,r)E'{r,p,q)/Eir,p,q)] 
Equation (|4]) can be integrated: 



As can be seen the Szekeres model is specified by 6 functions: 
M{r),k{r),tB{r),S{r),Q{r),P{r). However, by a choice of the coordinates, the 
number of independent functions can be reduced to 5. 

The equations of motion T"^;^ = are reduced to the continuity equation: 

p + pQ = 0, (7) 

where Q is the scalar of expansion and is equal to: 

(^(f ] = o ^(^'^) , <i>'(t,r)-$(^,r)a>^(t,r)/j>(t,r) 

^''^''^^ ^t,r)^^'{t,r)-^t,r)E'{r,p,q)/E{r,p,qy ^^ 

In the expanding Universe Q is positive so the density decreases. The structures 
which exist in the Universe, emerged either due to slower expansion of the space 
(formation of overdense regions) or due to faster expansion (formation of underdense 
regions). In the Friedmann limit R -^ ra, where a is the scale factor and Q —>■ 3Hq. 

The Szekeres model is known to have no symmetry.^ It is of great flexibility and 
wide application in cosmology''^ and in astrophysics,'*'^ and still it can be used as a 
model of many astronomical phenomena. In this paper it will be employed to study 
the evolution of cosmic structures in different environments. 
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2.1. Density contrast 

To compare the evolution of different models the change in their density contrast 
will be considered. Two different types of density contrast indicators are taken into 
account. The first one is the usual density contrast, S = p/ Pb — I7 the second one is 
the spatially invariant density contrast:^ 



5'iK — 



/i"^ dp dp ^ 



pi dx'-dxP '^^' ^^^ 

where / G M, and K E M\{0}. This family of the density contrast indicators can be 
considered as local or global depending on the size of E. Such a quantity not only 
describes the change of density but also the change of gradients and the volume of 
a perturbed region. So this density indicator describes the evolution of the whole 
region in a more sophisticated way than the S. Here only the case I = 2,K = 1/2 
will be considered. 

2.2. Model set-up 

To specify the model 5 functions of the radial coordinate need to be known. Let us 
define the radial coordinate as a value of $ at the initial instant to — 0.5 My after 
the big bang, i.e., r :— $(r, io)- 

Two of these functions will be tB{r) and M{r). Let us write the mass function 
in the following form: 

M{r) = Mo(r) + SM{r), (10) 

where Mq is the mass distribution as in the homogeneous universe, and SM is a mass 
correction, which can be either positive or negative. The SM is defined similarly as 
in the spherically symmetric case: 

6M{r)=4TT^ f drR{r,tofR'{r,to)dp{f), (11) 

c Jo 

where dp{r) is an arbitrary function chosen to specify the SM . Although Sp{r) is not 

the initial function of density fluctuations (since an initial density fluctuation is a 

function of all coordinates) it gives some estimation on the initial density fluctuation 

of the monopole density component. 

The next three functions are P{r), Q{r), S{r). All functions defining the model 
are presented in Table[TJ The numerical algorithm used to solve the Szekeres model's 
equations is presented in detail in Rcf. 11. 

The chosen background model is the homogeneous Friedmann model with the 
density: 

O 7j2 

Ph = nmX Per = 0.24 X -^. (12) 
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where the Hubble constant is Hq = 74 km s ^ Mpc ^. The cosmological constant, 
A, corresponds to flA = 0.76, where ft a = {l/3){c^A/H§). 

3. Models of a void and an adjacent galaxy supercluster 

In this section the evolution a void with an adjourning galaxy supercluster is in- 
vestigated. Although within the Szekeres model more than two structures can be 
described, such investigations of less complex cases are useful because they enable 
us to draw some general conclusions without going into too much detail, which could 
easily obscure the larger picture. 

3.1. Models with P' = = S", Q' ^ 

As mentioned above, ii P' = = S' = Q' the Szekeres model becomes the Lemaitre- 
Tolman model. Hence, the class of models considered in this subsection is the sim- 
plest generalisation of the spherically symmetric models. 

The double structure of a void and adjourning supercluster can be described in 
the Szekeres model in two different ways. The first alternative is when SM < 0, the 
second when SM > 0. Both these possibilities are examined here. 

3.1.1. Model specification 

The exact form of the functions used to define models 1 and 2 is presented in Table 
[H The density distributions of models 1 and 2 are presented in Fig. [TJ As can be seen 
the model with 5M < has the void in the center, and the adjacent supercluster 
has an elongated shape. It is the opposite in model 2. The overdense region at the 
origin is more compact than in model 1, and the adjacent element is the void. 

Table 1. The exact form of the functions used to define models 1-4. 
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-5A X exp[- 
1.14A X exp[ 
—5A X exp[- 
1.14A X exp[ 


-ie/8f] 

-(^/9)2] 

-{e/9f] 


1 
1 

_(0.A 
_£0.9 






0.55^0.4 


bln(l-|-^) xexp(-3A(?) 

cln(l + 0.2£) X exp(-3A^) 

0.33^°* 

0.33^°* 



Note: £ = r/Kpc, A = 10"^, b = -0.6, c = 1.45,, 



3.1.2. Evolution 



In this section we compare the evolution of the density contrast, S{t,r,0,(l)), and 
the 5'2,i/2(i, S) density indicator for models 1, 2, with the corresponding mod- 
els of a single void and the models of a single supercluster obtained within the 
Lemaitre-Tolman model. The Lemaitre-Tolman model is considered because within 
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Fig. 1. The present-day density distribution, p/pb- Fig. 1(a) presents model 1 {SAI < 0). Fig. 
1(b) presents model 2 {5M > 0). 



this model one can describe a single spherically symmetric structure. Such a com- 
parison can demonstrate how the evolution of a structure changes if there is another 
structure in its close neighbourhood. 

Fig. [21 presents the evolution of the density contrast of model 1 in comparison 
with the corresponding models obtained within the Lemaitre-Tolman model. The 
Lemaitre-Tolman model was specified by assuming the same condition as the ones in 
the Szekeres model at the initial instant. The local density contrast, S, is compared 
at the point of the maximal and minimal density value. Fig. [H^a) presents the 
evolution of the density contrast inside the void. As can be seen the bahaviour 
of the density contrast in both models is similar. This due to the conditions of 
regularity at the origin (for a detailed description of the regularity conditions at 
the origin see Ref. 4). These conditions imply that the origin in the quasispherical 
Szekeres and Lemaitre-Tolman model behaves like a Friedmann models Fig. [DJb) 
presents the evolution of the density contrast at the very center of the overdense 
region of the model 1 and the corresponding Lemaitre-Tolman model. The growth 
of density contrast in the Szekeres model is much faster than in the corresponding 
Lemaitre-Tolman model. The results of this comparison indicate that within the 
perturbed region of mass below the background mass {SM < 0) the evolution of 
underdensities does not change but the evolution of the overdense regions situated 
at the edge of the underdense regions is much faster than the similar evolution of 
isolated structures. 

The evolution of the density contrast of model 2 {SM > 0) is presented in Fig. 
[31 the evolution of the density contrast at the point of minimal density is depicted 
in Fig. El^a), and the evolution at the origin is depicted in Fig. ^h). Similarly as 
in model 1, the evolution at the origin in the Szekeres model and in the Lemaitre- 
Tolman model are very much alike. The evolution of the void, however, is slower 
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within the Szekeres model than it is in the Lemaitre-Tolman model. This implies 
that single, isolated voids evolve much faster than the ones which are in the neigh- 
borhood of large overdensities where the mass of the perturbed region is above the 
background mass {SAI > 0). 





time [10" y] 
(b) 

Fig. 2. The evolution of the density contrast inside the void (a), and inside the supercluster (b) 
for model 1 {SM < 0). The curve SZ presents the evolution within the Szekeres model; curve LT 
presents the evolution within the Lemaitre-Tolman model. 
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Fig. 3. The evolution of the density contrast inside the void (a), and inside the supercluster (b) 
for model 2 {SM > 0). The curve SZ presents the evolution within the Szekeres model; curve LT 
presents the evolution within the Lemaitre-Tolman model. 



3.2. Models with P' ^ ^ S", Q' ^ 

In this section models of non-constant P, Q and S are investigated. The evolution 
of these models is compared with the evolution of models which were considered in 
SecO 
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3.2.1. Models specification 

The exact form of the functions used to define models 3 and 4 is presented in 
Table [ij Fig. [4] presents the comparison of the present day density distribution in 
models 1 and 3 in colour coded diagrams. It presents the vertical cross-sections of 
the considered structures. Fig. ID^a) presents the vertical cross-section through the 
surface oi (j) = 7r/2 and Fig. [D^b) presents the cross section through the surface 
of « 7r/6. The comprehensive study of the vertical and horizontal cross-sections 
of similar models was presented in Ref. 12. Fig. [5] also presents the vertical cross- 
sections of models 2 and 4. As can be seen, both structures appear to be similar 
but, in comparison with model 1, in model 3 the additional component is moved 
down and right. Model 4 on the other hand presents the structure moved down and 
right in comparison with model 2. 
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Fig. 4. The present— day colour coded density distribution, p/ph, of models with 5M < 0. (a) 
presents the model with P' = S' = 0. (b) presents the model with P' j^ j^ S' . 



3.2.2. Evolution 

The evolutions of the density contrast inside the voids and superclusters of models 
3 and 1 are very similar which needn't be suprising as model 3 has the same 5{r) 
as model 1. Also, the evolutions of the corresponding density contrasts of model 4 
and 2 are similar. The functions S, P, Q were chosen so they reproduce the same 
shape of current structures and the same density contrast inside them — that is 
why the evolution of a local density contrast is comparable for models 1 and 3, 
and for models 2 and 4. However, it is not clear whether or not the evolution of 
52,1/2 is comparable too. When the functions S,P,Q are not constant, the axis of 
a density dipole changes. Also, the volume of the perturbed region as well as the 
density gradients can be different. So it may be interesting to compare the evolution 
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Fig. 5. The present— day colour coded density distribution, p/pbi of models with SM > 0. (a) 
model with P' = S' = 0. (b) model with P' ^Q^ S'. 



of the whole perturbed underdense and overdense regions of models 1, 2, 3, and 4. 
Fig. [6] presents the comparison of evolution of 52,i/2 for models 1-4. The primed 
letters denote models of S" 7^ 7^ P', Q' ^ 0. As can be seen the evolution of 82^1/2 
for all these models is also comparable. These results imply that the evolution in 
the quasispherical Szekeres model does not depend on the position of the dipole 
component. As long as the shape and density contrast of the analysed models are 
similar, such models evolve in a very similar way. 
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Fig. 6. Comparison of 52, 1/2 for models with SM < (a) and with 5M < (b). C corresponds 
to a "supercluster" — an overdense region, and V corresponds to a "void" — underdense region. 
Primes denote models with P' ^ ^ S' . Since the value of 5ik depends on units, the results 
presented here were normalised so they are now of order of unity. 
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4. The role of expansion 

The faster or slower evolution rate of the previously presented models is reflected 
by their current expansion rate. As has been shown above, the evolution does not 
depend on a relative position of the dipole component (evolution of models 1 and 3 
is similar). Thus, let us focus on model 1 and model 2 only. 

Fig. [7] presents the ratio, Ssz/0O; of the expansion parameter in the considered 
Szekeres models to the expansion parameter in the homogeneous background. As can 
be seen, model 1, with SM < 0, has a larger amplitude of this ratio, and the evolution 
of a supercluster in this model is much faster than in the corresponding Lemaitre- 
Tolman model. On the other hand, model 2 {SM > 0) has smaller amplitude of 
the Osz/00 ratio and within model 2 the evolution of the density contrast inside 
the void was much slower than in the Lemaitre-Tolman model. So clearly the rate 
of the evolution is connected with the rate of the expansion. This conclusion is 
also supported by the continuity equation [Eq. ([7|)]. Higher mass in the perturbed 
region slows down the expansion rate — this is a condition hindering the evolution 
of cosmic voids. On the other hand, if the mass of perturbed region is below the 
background mass, such region expands much faster than the background, leading 
to the formation of the large underdense regions. Such large voids enhance the 
formation of large elongated overdensities (walls) formed at the edges of voids. 
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Fig. 7. The ©sz/00 ratio, (a) presents the ratio of model 1. (b) presents the ratio of model 2. 



5. Conclusions 

The galaxy redshift surveys show that the Universe is patchy with various structures. 
These structures include small voids among compact superclusters and large voids 
surrounded by large walls or long filaments. 

The evolution of these cosmic structures in different environments in the quasi- 
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spherical Szckcres model was investigated. The Szekeres model is the most complex 
spatially inhomogeneous exact solutions of the Einstein field equations, and it is 
of great use in cosmology. Since it is an exact solution of Einstein's equations, it 
enables us to investigate the evolution of cosmic structures without such approxi- 
mations as linearity and small amplitude of density contrast. Moreover, the Szekeres 
model is flexible enough to describe more than one structure. 

Having investigated various models with two or three structures within one 
frame it may be concluded that the evolution of the cosmic structures depends 
on the environment. In perturbed a region whose mass is below the background 
mass the amplitude of the expansion's fluctuations is large and as can be seen from 
the continuity equation [Eq. ([7])] , such conditions enhance the evolution of cosmic 
structures. 

The analyses presented in this paper indicate that small voids among large 
overdense regions do not evolve as fast as the large voids do. This is because the ex- 
pansion of the space is faster inside large voids than inside smaller voids. Moreover, 
this higher expansion rate inside the large voids leads to the formation of large and 
elongated structures such as walls and filaments which emerge at the edges of these 
large voids. 
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